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Abstract. Motivated by recent work on the structure of the singularity in inhomogeneous 
Tolman-Bondi collapse models, we investigate the behaviour of null geodesics in the particular 
case where the collapse is self-similar. The presence of the homothetic Killing vector field implies 
that the geodesic equation can be described by an integrable Hamiltonian system, and exploiting 
this fact we provide a full qualitative picture for its phase flow. 


1. Introduction 

Our late colleague and friend Victor Villanueva had broad interests in theoretical physics, 
and although he was considered a particle physicist, he also made important contributions to the 
field of general relativity. In particular, Victor worked on the geodesic deviation equations 
for relativistic spinning particles and possible applications to the detection of gravitational 
radiation HI Ells], a topic that is tightly related to upcoming experiments regarding the detection 
of gravitational waves through the use of pulsar timing arrays j^. More recently, Victor also 
worked on the Hamiltonian formulation for higher-dimensional black holes |5|. We dedicate 
the present work to him. We are confident he would have shared our curiosity in the problem 
analyzed in this article. 

It is well known that a Tolman-Bondi spacetime, describing the collapse of a spherical dust 
cloud, admits shell-focusing singularities, a portion of which may be null and visible to local 
observers niziiHiEiiiniiiiiini. Furthermore, for suitable initial data, part of this null singularity 
may also be visible to observers in the asymptotic region, see for instance Bin], and in this case 
there is a Cauchy horizon which extends all the way to future null infinity. Although the Tolman- 
Bondi metrics are spherically symmetric and are known in closed explicit form, still establishing 
the existence of future-directed null geodesics emanating from the central singularity is a problem 
of considerable mathematical complexity. So far, sufficient conditions upon the initial data have 
been found that in turn guarantee the existence of radial null geodesics emanating from the 
central singularity B mil El]. 

The visibility of the null singularity to local or asymptotic observers has raised questions 
regarding the validity of the strong and weak cosmic censorship conjectures within Einstein’s 
general theory of relativity. It is not the intention of the present work to discuss this delicate 


and difficult matter, instead our intention is to provide a systematic discussion for the null 
geodesics (with and without angular momentum) emanating from or terminating at the central 
singularity. 

In order to achieve this goal, here we restrict ourselves to the family of self-similar Tolman- 
Bondi collapse models. For an introduction to such models, see for instance naiiaiisiiinKiTiiiH] 
as well as the review in |19| , and for questions related to the stability of the Cauchy horizon and 
cosmic censorship based on these models see [iiiEniEiiiasiEs]. Our motivation for restricting our 
attention to the self-similar case stems from the fact that the extra symmetry of the background 
geometry allows us to have a complete understanding on the behaviour of radial and non-radial 
null geodesics emanating from or terminating at the central singularity. This property is very 
welcome for a number of independent reasons. 

Some time ago, Mena and Nolan m have shown that any marginally bound nakedly 
singular Tolman-Bondi dust collapse model admits future-directed non-radial null geodesics 
emanating from the central singularity. At hrst sight, this result might appear unexpected 
and counterintuitive, since by conservation of angular momentum one expects that no causal 
geodesics possessing non-vanishing angular momenta can pass through the center of the cloud. 
As long as the center is regular, this expectation is indeed fulhlled. However, when the singularity 
forms, the singular nature of the center transcends the restrictions imposed by the law of angular 
momentum conservation, as Mena and Nolan show at least for the case of marginally bound 
collapse. In the present work, we determine the class of all possible (radial and non-radial) null 
geodesics emanating from the central singularity in the self-similar case, thus complementing 
the results in mills]. For an extension of these results to the more generic family of bounded 
nakedly singular Tolman-Bondi models, see our recent work |24) . 

There is another reason for which non-radial geodesics emanating from the central singularity 
play an important role, as we noticed in m In that work, via a combination of analytical and 
numerical computations, evidence is presented to support the view that these non-radial null 
geodesics emanating from the singularity shape the shadow that a collapsing cloud cast upon 
the image perceived by an asymptotic observer of an external source illuminating the collapsing 
cloud. Accordingly, these geodesics may herald to an observer in the asymptotic region the 
formation of a Cauchy horizon. 

The structure of the present paper is as follows. In the next section, we introduce the family of 
self-similar Tolman-Bondi spacetimes and discuss briefly their basic properties. In the following 
section, we discuss the important family of radially in- and outgoing null geodesics and explain 
the causal structure of the spacetime. In the subsequent section, we set up a Hamiltonian 
formalism describing the non-radial null geodesics and give a full qualitative description of the 
associated Hamiltonian flow. Finally, in the last section a summary of the results is presented 
and possible applications are discussed. 

2. Self-similar Tolman-Bondi spacetimes 

We begin this section reviewing the self-similar Tolman-Bondi collapse. The metric describing 
the collapsing spacetime is 

g =—dr'^F‘^{x)dR^R^S‘^{x) , x (1) 

R 

where r G M is the proper time along the flow lines of the dust particles, i? > 0 is a comoving 
radial coordinate labeling the dust shells, and (i?, ip) are standard angular coordinates on the 
two-sphere, while the functions F and S are dehned by 

-j _ ^ 


X < 1/A, 


( 2 ) 



with A a positive parameter. Denoting by m{R) the Misner-Sharp mass function, the parameter 
A characterizes the initial compactness ratio of each dust shell, 


2m{R) _ 4 2 
R 9^ 


which turns out to be independent of R in this model. The areal radius r(T, R) = RS{x) satishes 


1 /9r(r, i?)\^ m{R) 

2 \ (9r J r{T, R) 


( 3 ) 


implying that the collapsing shells have zero total energy. Furthermore, dr/dR = F(x) > 0 is 
positive, excluding the formation of shell-crossing singularities. 

The spacetime described by Eq. Q has a shell-focusing singularity at x = 1/A, where r/R 
vanishes and the density and the curvature blow up. Depending upon the range of the parameter 
A, this singularity may be visible to local observers (see Lemma below). In this work, we focus 
on the case where the shell-focusing singularity is visible and provide a detailed analysis for the 
behaviour of radial and non-radial null geodesics emanating from it. 

A particular property of the metric Q is the presence of the homothetic Killing vector held 




( 4 ) 


such that the Lie derivative of g along the integral curves of ^ scales according to = 2g. The 
existence of ^ yields a conserved quantity for the equations of motion describing null geodesics: 

Proposition 1. Let ^ be a homothetic Killing vector field admitted by a Lorentz manifold (M, g), 
and consider an affinely parametrized null geodesic 7 in (M, g) with tangent vector field p. Then, 
the quantity C := g(p,0 remains constant along 7. 


Proof. Using local coordinates (x^) on (M, g), we hnd 

The hrst term on the right-hand side vanishes because p is geodesic and affinely parametrized, 
while the second term vanishes because = £^g^v = “^g^u and p is null. □ 

Notice that the norm of ^ is 

g(C,0 = -R^{x^ - F‘^{x)), 

implying that ^ is timelike in the region where x‘^ — F‘^{x) > 0 but spacelike where x^ — F‘^{x) < 0. 
Therefore, C cannot be interpreted as a conserved “energy" in an obvious way. Nevertheless, 
as we will see, the existence of the conserved quantity C greatly facilitates the analysis of the 
geodesic flow with non-vanishing angular momentum. 


3. Radial null geodesics 

In this section, we analyze the behaviour of radial null geodesics. From Eq. Q, it follows 
that these geodesics are determined by 

dr = ±F{x)dR, 

where the -l- (—) sign refers to outgoing (ingoing) radial null geodesics. In terms of the variables 
X = t/R and s := — log(i?) this equation takes the form 

dx 

— = Kt(x) := X =F F(x), x < 1/A. (5) 

ds 

The qualitative behaviour of the solutions of these equations can be understood from the 
properties of the functions Y±, which are summarized in the following lemma. 






Figure 1. A plot of the two functions 1+ (blue-continuous) and IL (red-dashed) for the 
parameter value A = 0.5. In this case, Jq — —1.0198, J\ ~ 1.0588 and J 2 ~ 1.9088. 


Lemma 1. The functions Y± : (— 00 ,1/A) —>■ M defined in Eq. 0 satisfy the following properties 
(cf. Fig. 

(a) lim Kt(x) = — 00 , lim Kt(x) = =Foo. 

x^-00 x^l/X 

(b) Y- has a single root Jq satisfying Jq < 0. 

(c) Regarding the function Y^, let 

X* := ^(V3 - l)(3\/3 - 5)^/3 ~ 0.638. 

For A > A*, 1+ is strictly negative, for A = A*, 1+ has a single root at Ji = (\/3 — l)(3-v/3 — 
5)“^/^, while for 0 < A < A*, 1+ has precisely two roots Ji, J 2 satisfying 0 < Ji < J 2 < 1/A. 


Proof, (a) follows by inspection. For (b) and (c) the following expressions are helpful: 

F'{x) = ^^^3;(1 - F"{x) = ^A^ ^1 -h (1 - Ax)“’^/^. 


In order to prove (b), we hrst note that the existence of the root follows from the intermediate 
value theorem. There are no positive roots because F_ (0) = 1 and Y- is increasing in the interval 
(0,1/A). If Jq < 0 is a negative zero of YL, then it is not difficult to show that Y!_{Jq) > 0, 
implying the uniqueness of Jq. 

Finally, in order to prove (c), we note that F+(0) = —1, ^/.(O) = 1, and that 1+ is increasing 
for X < 0 and concave for x G [0,1/A). Therefore, 1+ has a global maximum at some point 
Jmax £ (0,1/A) and depending on whether the value of y+ at this maximum is negative or not, 
the function Y/_ has no roots, or it has precisely two roots J 2 > Ji in the interval (0,1/A). If Ji 
is a root of y+ then, zi := AJi G (0,1) satishes zi = A(1 — 2;i/3)(l — zi)~^/‘^, and 


A = P{zi) 


Zl(l-Zi)V3 



The positive function P : (0,1) —>■ M converges to zero for 2 : —)> 0 or z —>■ 1 and it has a maximum 
at zi = Zi = 3(2 — \/3), where 

A* = P{zl) = 1{V3- l)(3\/3 - 5)^/3^ 









Therefore, the function 1+ has no roots when A > A*, it has a degenerate root at Ji = J 2 = 
z\/\* = (\/3 — l)(3-v/3 — 5)“^/^ when A = A*, and it has two roots when 0 < A < A*. □ 

For the following, we assume that the parameter A in Eq. Q has been chosen so that it lies 
below the critical value A* i.e. 0 < A < A*, and thus the function 1+ has precisely two zeros 
denoted by Ji < J 2 . Consequently, in the outgoing case, the system ([^ has two critical points at 
X = Ji and X = J 2 . It follows from the behaviour of 1+ that x = Ji is unstable, while x = J 2 is 
an attractor: any solution with initial data x(0) G (Ji, 1/A) converges to J 2 as s —)> 00 (-R —)• 0), 
while any solution with x(0) G (— 00 , Ji) converges to —00 when s —)• 00 (i? —)• 0). The solution 
is implicitly determined by the equation 


a:(s) 

s-so= [ -(6) 

J X- F{x) 
xlo) 

It follows that any outgoing radial null geodesic that passes through a point (r, R) such that 
Ji < t/R < 1/A, R > 0, emanates from the central singularity {t,R) = (0,0), see Fig. On 
the other hand, the outgoing radial null geodesics passing through a point (r, R) with t/R < Ji 
and R > 0 emanate from a point (tq < 0, 0) on the regular center. 




Figure 2. Spacetime diagrams illustrating the causal structure of self-similar dust collapse. 
Left panel: The location of the singularity and the particular radial null rays x = Jq , Ji, J 2 in the 
T — R plane. The green solid lines with the arrows represent future-directed radial null geodesics 
emanating from the central singularity. Right panel: The location of the same objects in the 
conformal diagram. 

Regarding the case of the ingoing radial null geodesics, it follows from Eq. ([^ and the 
behaviour of IT, that x = Jq is an unstable critical point. This implies that any ingoing radial 
null geodesic passing through a point (r, R) with t/R < Jq and ii > 0 terminates at the regular 
center, while any ingoing radial null geodesic passing through a point (t, i?) with t/R> Jq and 
R > 0 terminates at the shell-focussing singularity with finite value of R, as can be seen from 
the convergence of the integral 

x{s) 

f dx 

J X + F{x )' 

x(0) 


( 7 ) 










when x(0) > Jq and x(s) —)• 1/A. 

Therefore, there is a unique ingoing radial null geodesic terminating at the central singularity 
(r, R) = (0, 0), while there are infinitely many outgoing radial null geodesics emanating from it. 
The earliest one, x = R, describes the Cauchy horizon, see Fig. 


4. Hamiltonian description of the null geodesic flow 

In the previous section, we reviewed the behaviour of radial in- and outgoing null geodesics. In 
this section, we shall appeal to the powerful Hamiltonian formalism to treat the more general case 
of non-radial null geodesics in the background metric described in Eq. Q. Via this formalism, 
based on the rotational symmetry and the presence of the homothetic Killing vector field, we 
show that the description of the complete set of null geodesics is reduced to an effective one¬ 
dimensional Hamiltonian system. 

We recall that geodesic motion is described by the Hamiltonian 


and the null geodesics are those trajectories along which H = 0. Because of spherical symmetry, 
it is sufficient to restrict ourselves to the equatorial plane = 7r/2 and moreover it is convenient 
to express H{x,p) in terms of the local coordinates (x, s, i?, p) where x = t/R and s = — \og{R). 
Relative to these coordinates, the homothetic Killing vector held assumes the simple form 



and the inverse metric on the equatorial plane reads 

d d 1/(9 d \ / d d \ Id d ' 

dxdx F'^{x) \ds ^dx) \ds ^dx) S'^{x) dp^ dp 

Since conformal transformations do not affect the null geodesics as trajectories, hereafter, we 
discard the conformal factor The resulting Hamiltonian takes the form 



H{s,x,p,ps,ps,p^) 


1 

2 


2 {Ps+ XPxf P% 
F 2 (x) 52 


The variables s and p are cyclic, and thus the following two quantities 


C :=-ps, i-=P,p, 


( 8 ) 


are preserved along the flow. The quantity C is the constant of motion associated with the 
homothetic Killing vector held see Proposition!^ while i is the constant of motion associated 
with the rotational Killing vector held djdp. As a consequence of these remarks, the Hamiltonian 
can be reduced to the effective one-dimensional system 


Hc,e.{x,Px) = ^ -pI + 


{xpx - C)2 ^ /2 


F2(x) ' 52 (x) ■ 

This effective Hamiltonian describes the motion in the x — px plane: 
dHci 


X = 


Px = 


dp, 

dHc,i 

dx 


= -Px + 


F2(x 


-{xpx - C), 


(9) 


( 10 ) 


Px . , {xPx - Cf , f , 

[xpx -C) -\ -—7WW —F (x) -F „ S (x). 


F2(x) 


F3(x) 


53 (x) 


( 11 ) 




















and we are only interested in those trajectories for which Hc^i = 0. Eliminating px in Eq. (10) 
by using Hc^i = 0 , we obtain 


+ Vc,^{x) = 0, Vc,,{x) := " C^S\x)\ , (12) 

which is the equation for a point particle with zero energy moving in one dimension. 

Once the motion in the x — Px plane has been determined, the additional degrees of freedom 
describing the flow can be obtained by integrating the equations 

s 


dH _ xpx — C 
dps F‘^{x) ’ 
dH _ i 
dp^ S‘^{x)' 


(13) 

(14) 


4-1. The radial case (£ = 0) 

As a first application of the formalism developed in this section, here we briefly reconsider 
the case of radial null geodesics. The equations for these geodesics are obtained by setting £ = 0 
in Eq. (12) which yields x = ±C/F{x). On the other hand, the constraint He/ = 0 yields 


[(x - F{x))px - C] [(x + F{x))px - C] = 0 . 
Assume first that T^(x) 7 ^ x^. Then, we obtain from this the solutions 

C 

Px — 

X =F F{x) 

Introducing this relation into Eqs. ( |10|13[ ) we obtain 

C .Cl 


X = ± 


F(x)’ 


s = ± 


F{x) X =F F{x)' 


(15) 


If C* = 0, then x = s = 0 and Px = Ps = Pip = 0, so we obtain a trivial geodesic. If C 7 ^ 0, we 
can divide x by s and we recover Eq. 

When T^(x) = x^, we obtain from Hq/ = 0 either C = 0 or = Cj{2x). In the first case, 
X = 0 and s = Pxjx, and thus we recover the critical radial null geodesics x = Ja, a = 0 , 1,2 
discussed in the previous section. In the second case, we obtain x = —Cjx and s = —C'/(2x^), 
which is consistent with the limit of Eq. (15) when F{x) — >• ^x. 


4 . 2 . The non-radial case with (7 = 0 

When (7 = 0 and ^ 7 ^ 0, it follows from Eq. (12) that the motion is restricted to the region 
where x^ — T^(x) < 0, implying that ^ is spacelike or null, and in this case, the Hamiltonian 
constraint He/ = 0 simplifies to 


{x^-F\x))pl + i 


•FM 

52(x) 


= 0 . 


Since CF^jS^ is positive, neither x^ — F^{x) nor px can vanish, and it follows that ^ is spacelike 
and that the motion is confined to the regions Jq < x < Ji or J 2 < x < 1/A. Dividing Eq. (13) 
by Eq. (10) we obtain 

ds X 

dx F2(x) — x2 ‘ 

















Integration leads to 


s = So - 


x{s) 


xdx 


F‘^{x) — ’ 


XO 


(16) 


which determines the trajectory of the light ray in the x — s plane, for given initial data (xq, sq) 
with So G M and xq lying either in the interval (Jo, Ji) or in the interval (J 2 ,1/A). In the first 
case, s —)> —00 {R —)> 00) as x(s) —)• Jq or x(s) —>■ Ji, and thus we obtain a null geodesic which 
asymptotes to the Cauchy horizon x = Ji in the future and to x = Jo in the past. In the second 
case, s — )> 00 (i? — )• 0) as x(s) —)> J 2 , while s and R converge to a finite value when x(s) —)• 1/A, 
so in this case we have a null geodesic emanating from the central singularity which terminates 
at the spacelike portion of the singularity. 


4-3. The generic case (i C 7 ^ 0/ 

Finally, we treat the generic case of non-radial null geodesics with constants of motion C and 
I different from zero. In this case, the motion in the x-direction is restricted to the set of points 
for which (see Eq. (dl)) 


W{x) : = 


g(C,6 _ 

J.2 


where we have introduced the “impact parameter" /3 : 
in the following lemma. 


-FHx) 1 
52 (x) -, 02 ’ 


(17) 


= ijC. The structure of this set is discussed 


Lemma 2. Consider the effective potential W : (— 00 ,1/A) —>■ M defined by Eq. & For each 
(3 0 let Ig be the set consisting of those points x for which W{x) < 1/fifi. Then, there exists 

fic > 0 depending on A such that 


Id 


(xq, 1 /A) for some xq < Jo, < /?/, 

(xo, 1/A) \ {xi} for some xq < Jq and Ji < xi < J 2 , f3‘^ = /3/, 

(xo,xi) U (x 2 ,1/A) for some xq < Jo and Ji < xi < X 2 < J 2 , > fic 


(18) 


The allowed region for the motion in the x-direction is the closure, Ig in (—00, 1/A), of the set 

Id- 

Proof. The functiou W : (—00, 1/A) —M satisfies lima;-,.-oo W{x) = +00 and lim 2 ,_,.i/;, W{x) = 
—00, and as a consequence of Lemma[^it is positive on the intervals (—00, Jq) and (Ji, J 2 ) and 
negative on the intervals (Jq, Ji) and (J 2 ,l/A), with simple roots at Jq, Ji, J 2 . Furthermore, a 
straightforward calculation reveals that 


W'{x) = 2 ^1 — (1 — Ax) X — Ax) 


(19) 


implying that W is decreasing for x < 0 . Hence, for any > 0 there exists a unique xq < Jo 
such that l//3^ = ly(xo). Therefore, Ig C (xo,l/A). If is small enough, it is clear that the 
inequality 0 is satisfied for all x G (xo,l/A), because W{x) —>■ —00 as x —)• 1/A, so in this 
case Ip = (xo, 1/A). On the other hand, for large /3^, the inequality Eq. ( |l7[ ) cannot hold for all 
X G (Ji, J 2 ) since W is positive on this interval. 

We now claim that the function W' has a unique root Xc on the interval (Ji, J 2 ), corresponding 
to a maximum of W on (Ji, Jfi). The statement of the lemma then follows with fic ■= W (xc)~^^^. 
In order to prove the claim, we consider the function U : (Ji, J 2 ) —)• M, x 1 —x —(2A/3)(1 —Ax)“^/^ 
which determines the sign of W'{x), see Eq. (19). From the behaviour of the functiou x^ — F‘^{x) 








that follows from Lemma we know that W'{Ji) > 0 and W'{J 2 ) < 0, implying that U is 
positive near x = Ji and negative near x = J 2 . By the intermediate valne theorem, there exists 
a point xi G (Ji, J 2 ) snch that U{x) > 0 for all x G (Ji,xi] and U'{xi) < 0. Since U is concave, 
it follows that U'{x) < 0 for all x G (xi, J 2 ) which implies that U has a nniqne zero on (Ji, J 2 ). 
This conclndes the proof of the lemma. □ 

Restricting onrselves to the allowed region Ip, we nse = 0 to obtain 


Px = 


X 

C_ 

2x 


[x ± F{x)Q 0 {x)] , F2(x)/x2, 
F2(x)=x2 


( 20 ) 


with 


Qpix) := Vl-= y 1 - 

In the hrst case, when F‘^{x) 7 ^ x^, we obtain from Eqs. ( 10|l3 ), 

C ^ . . . C xQp{x)±F{x) 




In the second case. 


C 

X = -, 

X 


s = ± 


C 


F{x) x^ —F2 (x) 




52 (x) 


( 21 ) 


( 22 ) 


which is seen to be the limit of Eq. (21) when F{x) —)■ q^x. We also see that Eq. (21) rednces to 
the corresponding eqnation (15) in the radial case when /3 —?■ 0. 

For the analysis below, the following identity is worth noticing: 

{xQp{x)f - F^{x) = (x^ - F^{x)) (^1 - 

Using this identity and Eqs. (2lp2) we obtain 


-r = ^0±\x) '■= 
ds 1 - / 32 ^ 


[xQp{x)^ Fix)) = Qp{x) 


S^ix) 


xQp{x) ± F{x) ’ 


(23) 


The fnnctions Yp^± determine the trajectories of non-radial nnll geodesics in the t — R plane. 
Note that these fnnctions rednce to the fnnctions Y± describing radial nnll geodesics when /3 —)■ 0, 
see Eq. (|^. The qnalitative properties of the fnnctions Yp^± are snmmarized in the next lemma 
and illnstrated in Fig. 

Lemma 3. Denote hy Di < i) < D 2 the two roots of the equation (I'^x'^ = S‘^{x). Then, Yp^j,. in 
Eq. (23) yields a well-defined smooth function Yp^^ : Ip \ {Hi} —)• M. It has a first-order pole at 
X = Di, is positive on (xq, Hi) U (Ji, J 2 ) and negative on (Hi, Ji) U (J 2 ,1/A). 

Similarly, Yp_ in Eq. 1^2!^ yields a well-defined smooth function Yp_ : Ip \ {H 2 } —)• M which 
has a first-order pole at x = D 2 , is negative on (xq, Jq) U (H 2 ,1/A) and positive on (Jq, H 2 ). 

Remark: Note that Hi and H 2 lie inside the region Ip, since 


1 

w 


F^Dg) ^ 1 
S^Da) / 32 ’ 


a = 1 , 2 . 


W{Da) 























Figure 3. A plot of the two functions 5^,+ and Yjs- for the parameter value A = 0.5. Left 
panel: /3 = 0.2, which lies below the critical value /3c. The two graphs connect at the point xq 
which lies around x = —35.2. Right panel: /3 = 0.9, which lies above the critical value [3c- In 
this case, it is clearly visible how the graphs of the two functions To. 9 ,± connect to each other at 
the turning points x = xq, x = xi and x = X2- 


Proof. The proof of Lemma follows directly from the two representations of in Eq. (23) 
and the known behaviour of the function — F‘^{x). 


□ 


From the qualitative behaviour of the functions described in Lemma one finds the 
following properties of the null geodesics: consider first a solution of dx/ds = with 

initial data at (x, s) = (xq, sq), where xq is the turning point defined in Lemma This solution 
extends through the pole at x = Di to x = Ji- Since 


x(s) 



xo 


the pole of T^,+ simply describes a turning point of s (or R). Thus, this solution describes a 
null geodesic emanating from the point (xq, sq) which asymptotes to the Cauchy horizon in the 
future. To describe the past of this geodesic, one has to resort to the function since xq is a 
turning point for the motion in the x-direction. From the properties of Y^_, one sees that the 
null geodesic asymptotes to the light ray x = Jq in the past. 

Next, consider a light ray described by Ts,+ passing through a point (x, s) with Ji < x < 1/A. 
For (3‘^ < /?/ the picture looks qualitatively the same as in the radial case. In particular, all such 
null geodesics emanate from the central singularity. Regarding the null geodesics described by 
Ya^_ passing through a point (x,s) with Jq < x < 1/A, they asymptote to x = Jq in the past, 
and terminate at the spacelike portion of the singularity, see Fig. 

Therefore, when /3^ < /?/, we have the following properties: all the null geodesics passing 
through a point (x, s) with xq < x < Jq bounce off at x = xq and asymptote to the Cauchy 
horizon in the future, while null geodesics passing through a point (x, s) with Jq < x < Ji 
terminate at the spacelike portion of the singularity. Finally, null geodesics passing through a 
point (x, s) with Ji < x < 1/A emanated from the central singularity (r, R) = (0, 0). 
















For as in the previous case, there are null geodesics which asymptote to x = Jq in 

the past direction and to x = Ji in the future direction, bouncing off at x = xq- However, the 
qualitative features of the remaining null geodesics, namely those that penetrate or lie to the 
future of the Cauchy horizon are different than in the previous case: the effect of having a high 
angular momentum such that is the appearance of the “forbidden" region xi < x < X 2 

(see Lemma 1^, which prevents the null geodesics originating from the far past r —)> —oo to reach 
the singularity. As increases, the gap between x = xi and x = X 2 increases, until xi —)• Ji 
and X 2 —>■ ^ 2 ) see Eq. ( 0 . 




Figure 4. t — R plane showing future-directed non-radial null light rays (green lines with 
arrows) propagating in the spacetime describing self-similar dust collapse. Left panel: < /3^. 

Right panel: /3^ > /3^. 


Finally, we comment on the behaviour of the azimuthal angle (p for the null geodesics 
emanating from the central singularity. In the cases observed so far, where 0 < /3^ 7 ^ ,0^, 
these geodesics have x —)• J 2 as they approach the singularity. It follows from Eqs. (14pT) that 


^ , /3 F{x) 

dx S^{x)Qi3{x)^ 


implying that ip has a hnite limit as x —?• J 2 . However, notice that in the particular case where 
/3 = ±/3c there exist null geodesics conhned to the timelike surface x = Xc = const., the value 
Xc corresponding the the local maximum of the effective potential W. In this particular case, it 
follows from Eqs. (|14|21|) that 

dp^ l 
ds 


implying that |</ 5 | —>■ 00 as the central singularity is approached. This situation is somehow 
analogue to the existence of unstable circular null geodesics in the Schwarzschild spacetime. 


5. Conclusions 

In this work, by appealing to the homothetic and spherical symmetry of a nakedly singular, 
self-similar Tolman-Bondi spacetime, we have presented a complete analysis of the behaviour 












of null geodesics on such spacetime. The results of this work confirm the existence of future- 
directed non-radial null geodesics emanating from the central singularity and are in agreement 
with those obtained by Mena and Nolan m and also with our more general results in |24) . 
The Hamiltonian treatment employed in this work has proven to be a very efficient technique to 
obtain many insights into the behaviour of null geodesics. For instance, through this technique, 
it became clear that the magnitude of the homothetic generator and the area of the two spheres 
define an effective potential W{x) and a critical impact parameter which plays a decisive 
role in the behaviour of radial and non-radial null geodesics as Lemmata and show. Besides 
the fact that they shed light on the structure and properties of the central singularity of a 
self-similar Tolman-Bondi spacetime, the results of the present work are important in another 
respect. They offer the possibility of studying the shadow that a naked singularity casts into the 
eyes of an asymptotic observer. This analysis would require the interior self-similar collapsing 
spacetime to be matched to an exterior Schwarzschild vacuum spacetime, a task that can always 
be accomplished satisfying standard junction conditions (see for instance |25] and references 
therein). Details of the modeling of the shadow of a naked singularity are discussed elsewhere |24) . 
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